Steady-state entanglement in a double-well Bose-Einstein condensate through 

coupling to a superconducting resonator 



O 



H. T. Ngi and Shih-I Chu^'^ 
^ Center for Quantum Science and Engineering, Department of Physics, 
National Taiwan University, Taipei 10617, Taiwan and 
'^Department of Chemistry, University of Kansas, Lawrence, Kansas 6604-5, USA 

(Dated: January 26, 2013) 

We consider a two-component Bose-Einstein condensate in a double-well potential, where the 
atoms are magnetically coupled to a single-mode of the microwave field inside a superconducting 
resonator. We find that the system has the different dark-state subspaces in the strong- and weak- 
tunneling regimes, respectively. In the limit of weak tunnel coupling, steady-state entanglement 
between the two spatially separated condensates can be generated by evolving to a mixture of dark 
states via the dissipation of the photon field. We show that the entanglement can be faithfully indi- 
cated by an entanglement witness. Long-lived entangled states are useful for quantum information 
processing with atom-chip devices. 
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I. INTRODUCTION 

Recently, the realization of a Bose-Einstein condensate 
(BEC) strongly coupled to the quantized photon field in 
an optical cavity has been shown [l], Q . This paves the 
way to study the interplay of atomic interactions and 
atom-photon interactions. For example, the novel quan- 
tum phase transition of a condensate coupled to a cav- 
ity has been demonstrated [Sj] • Strong atom-photon cou- 
pling is useful for quantum communications [41 such as 
the light-matter interface [l], [2| • 

Alternatively, strong coupling of ultracold atoms to a 
superconducting resonator has been recently proposed 
0. The two long-lived hyperfine states |e) = \F = 
2,TOF = 1) and \g) = \F = l,mF ^ -1) of *^Rb H 
are considered to be magnetically coupled to the mi- 
crowave field via their magnetic dipoles d, 0] ■ Since the 
high-Q superconducting resonator can be fabricated to a 
small mode volume ^ and the coupling strength can be 
greatly increased due to the collective enhancement 
The strong coupling of ultracold atoms in the microwave 
regime can be achieved Q . 

In this paper, we study a two-component BEC in a 
double- well potential , where all atoms are equally cou- 
pled to a single-mode of the microwave field inside a su- 
perconducting resonator. Two weakly linked condensates 
can be created in a magnetic double- well potential on an 
atom-chip [13, [HI or in an optical double-well poten- 
tial [11] • In fact, the tunneling dynamics between the 
atoms in two wells has been recently observed [l3l - [l5| . 
A double-well BEC coupled to an optical cavity has also 
been discussed in the literature [l6l - [l8| . However, the 
spontaneous emission rate of excited states used for opti- 
cal transitions in experiments ^i, i2] is much higher than 
the tunneling rate of the atoms between the two wells [Tsl - 
[T5j . Here we consider the two hyperfine states |e) and \q) 
of *^Rb with the transition frequency 2n x 6.8 GHz 6]. 
The coherence times [lll,!!^] of these hyperfine spin states 
(|e) and \g)) are much longer than both the timescales of 



tunneling and atom-photon interactions. Therefore, this 
system offers possibilities for the study of how the tunnel 
couplings between the two spatially separated conden- 
sates affect the atom-photon dynamics. 

We focus our investigation on the the system in the 
limits of the strong and weak tunnel couplings, respec- 
tively. We find that the system has the different dark- 
state subspaces [2l[ in these two tunneling regimes. In 
the weak-tunneling regime, the system has a family of 
dark states which can be used for producing quantum 
entanglement between the condensates. Here we propose 
to efficiently generate steady-state entanglement between 
the two spatially separated condensates by evolving to a 
mixture of dark states through the dissipation of the pho- 
ton field [2342^ . Note that our scheme does not require 
any adjustment of the tunneling strength. It is different 
to other methods ^] which depend on the strength of tun- 
nel couplings to generate entanglement. In addition, the 
entanglement generated between the two condensates can 
be used for the implementation of quantum state trans- 
fer 25|. This may be useful for quantum information 
processing with atom-chip devices [20| . 



Two-component BEC 




FIG. 1. (Color online) Schematic of a two-component BEC 
coupled to a single-mode of the photon field inside a supercon- 
ducting resonator. A two-component condensate is trapped 
in a double- well potential, and it is placed close to the surface 
of the superconducting resonator. The atoms are coupled to 
the magnetic field via their magnetic dipoles. The parame- 
ters L and w are the length and width of the superconducting 
resonator, respectively. 
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This paper is organized as follows: In Sec. II, we in- 
troduce the system of a two-component condensate in a 
double- well potential, and the two-level atoms are cou- 
pled to a superconducting resonator. In Sec. Ill, we de- 
rive the two effective Hamiltonians in the strong- and 
weak-tunneling regimes, respectively. In Sec. IV, we in- 
vestigate the dark-state subspaces and the atom-photon 
dynamics in the two tunneling limits. In Sec. V, we pro- 
vide a method to produce the steady-state entanglement 
between the two condensates in a double well. A sum- 
mary is given in Sec. VI. In Appendix A, we discuss the 
validity of the effective Hamiltonian in the strong tun- 
neling regime. 



II. SYSTEM 

We consider a two-component BEC being trapped in 
a double- well potential [9] , and the condensate is placed 
near the surface of a superconducting resonator as shown 
in Fig.[T] The atoms, with two internal states |e) and \g), 
are coupled to a single mode of the photon field via their 
magnetic dipoles. 



A. A two-component condensate trapped in a 
double-well potential 

We first introduce the system of a two-component con- 
densate in a one-dimensional (ID) double-well potential 
which can be described by the Hamiltonian as 



2ma dx"^ 

x-^o.{x)\'^o.{x) + 2U,gjdx-^\{x)'^l{x)-^g{x)-^,{x), 

(1) 

where ^'q. {x) is the field operator of the atoms for the in- 
ternal state I a) at the position x, and the indices a = g^e 
represent the ground and the excited states, respectively. 
Here ttIq is the mass of the atom in the state \a) and 
VvtYjjx) is the ID double- well potential which is given by 
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the effective ID interaction strengths between the inter-, 
and the intra-component condensates, as [1^ [27j 
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where C « 1.4603. The parameters Oq and aeg are the 
three-dimensional s-wave scattering lengths for the inter-, 
and the intra-component condensates. 

We adopt the two-mode approximation '28] such that 
the field operator (x) can be expanded in terms of the 
two localized mode functions and Uaj^ {x) as, 

'^a{x) = aLUabix) + ORUaRix), (5) 

where and an are the annihilator operators of the 
atoms in the state a = e, g for the left and right modes of 
the double-well potential, respectively. The Hamiltonian 
of the system [9|, within the two- mode approximation, 
can be written as 



^0 



hEeielcL + e^'^en) + hEg{g]^gR + glgL) 
-hJeielcR + e^CL) - hJg{g[gR + g\^gL) 
+hUee[ieleLf + {elcRf] + hUgg[{gigLf 
+{9R9Bf] + 2hUeg{e[eLglgL + e^ReRg^j^gR), 
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and j = L,R. The positive parameters Ea and Ja [29[ 
are the ground-state frequencies of the localized mode 
ctL,R, and the tunneling strengths between the two wells 
for the atoms in the states a. Here Ua and Uap are the 
two positive parameters which describe the inter- and 
intra-component interaction strengths, respectively. 



Vuwix) = Vd 



\Xo 



(2) 



where Vd is the barrier height and xq is the distance be- 
tween the two separate potential wells. The atoms are 
transversely confined in the y- and z-directions with the 
trap frequencies lu±. The size of the ground-state wave 
function in the transverse motion is a± — y/h/maUJ± 
[26l [27| , where me and nig are nearly equal. Since the 
transverse frequencies are much larger than the trap fre- 
quency in the x-direction, the transverse motions of the 
atoms are frozen out. The parameters Ua and Ueg are 



B. Atoms coupled to the photon field in a 
microwave cavity 

We consider that the atoms are coupled to a single- 
mode of the photon field via their magnetic dipoles [5j. 
Within the two-mode approximation, the Hamiltonian, 
describing the system of cavity field, the atoms and their 
interactions, is given by 

Hj = hwaO^a + nwo{e\eL + e^e_R) -f hg[a{e\gL + e^^fl) 
+H.C.], (11) 
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where Wa and a are the frequency and the annihilator 
operator of the single-mode of the photon field, and Wq is 
the transition frequency of the two internal states. Here 
we have assumed that the wavelength of the microwave 
field 1 cm) is much larger than the size of the conden- 
sate 10 iMm) [13, Eli- Therefore, all atoms are cou- 
pled to the photon field with the same coupling strength 
g = ^B\/^M)'^Ja/^^KV [7j], where /is is the Bohr magneton, 
/io is the vacuum permeability and V is the volume of the 
superconducting resonator. The coupling strength g can 
attain 1 kHz [3| if the volume V of the superconducting 
resonator is taken as Lxwxth ^ 1 cm x 10 /im x 200 nm 
05 111 J where L is the length, w is the width, and th is the 
thickness of the superconducting resonator. 



III. EFFECTIVE HAMILTONIANS IN STRONG 
AND WEAK TUNNELING REGIMES: LOW 
ATOMIC EXCITATIONS 

We will derive the effective Hamiltonians of the system 
in the limits of strong and weak tunnel couplings, respec- 
tively, where a few atomic excitations are only involved. 
Let us first write the total Hamiltonian of the system as 

H = hhJaO^a + hhjQ{e\eL + e^i^en) - hJe{e\eR + g^cl) 
-fiJgialgB. + 9r9l) + hUee[{e\eLf + {e^Renf] 

+eflefl5fl5fl) + ^9[a{e\9L + e\9R) + H.c] . (12) 

The total number of atoms N is conserved. We have 
omitted the constant term Et^N for a symmetric double 
well, where Ea « Eq for the two masses rUf, and rUg being 
equal. It is convenient to work in the rotating frame by 
applying the unitary transformation to the Hamiltonian 
H in Eq. (|12l) . where the unitary operator U{t) is 

U{t) = cxp [—iuja{a^a + e\eL + e\eR)t]. (13) 

The transformed Hamiltonian becomes 

H' = hA{e\eL + e^Cfl) - hJe{e\eR -f ej^ei) - hJg{g\gR 
+9r9l) + hUee[{eWf + (4efl)'] + f^Ugg[{g{gLf 
Hg^^Rgnf] + '2KUeg{e\eLglgL + e'^RCRgl^gR) 
+hg[a{e[gL + ej^gfl) + H.c] (14) 

where A = luq — uja '^s the detuning between the frequen- 
cies of the photon field and the two internal states. 

In the strong tunneling regime, the tunnel coupling 
is dominant and the strength of atom- atom interactions 
is relatively weak. On the contrary, in the weak tunnel- 
ing regime, the atom-atom interactions become dominant 
and the tunneling strength is negligible. We will show 
that these two cases exhibit the different behaviours in 
the atom-photon dynamics. We will provide derivations 
of the two effective Hamiltonians in the two tunneling 
limits in the following subsections. 



A. Strong-tunneling regime 

In the limit of the strong tunnel coupling, the tunneling 
strengths are much larger than the strengths of the atom- 
atom interactions, i.e., Je, Jg 3> Ue,Ug,Ueg- The total 
Hamiltonian of the system can be approximated as 

Hi = hA{e\eL + e\eR) ~ hJe{e\eR -I- e^ei,) 

-f^Jg{9i9R + 9r9l) + hg[a{e[gL + e^RgR) + H.c.]. 

(15) 

We have neglected the terms of the atom-atom interac- 
tions in this Hamiltonian. 

The symmetric and asymmetric modes g± and e± can 
be related to the localized modes as 

5± = -^(5L ±5i?,), (16) 

e± = -^(eiieij). (17) 

The Hamiltonian is then transformed as 

H[ = h{A - Je)e^e+ + h{A + Je)eLe_ - hjg{glg+ 

-9-9-) + f^9{ae\_g+ + H.c.) -I- hg{ae^_g_ + H.c). 

(18) 

Here the atoms are in symmetric (asymmetric) mode if 
they are populated in the states g'ji 1 0) + or e*^ 1 0) + (.g*i 1 0) _ 
or e'l|0)_), where |0)+ (|0)-) is the vacuum state of the 
symmetric (asymmetric) mode and fc is a non-negative 
integer. 

We consider the system to be initially prepared in the 
ground state in the limit of strong tunnel coupling, i.e., 
the ground state of the symmetric mode. The ground 
state can be obtained by applying the operator (g^)^ to 
the vacuum state |0)+ of the symmetric mode, i.e., 

1*1(0)) ^^(g;)^|0) + , (19) 

where N is the total number of atoms. Note that the 
atoms in the symmetric and asymmetric modes are inde- 
pendently coupled to the photon field in Eq. ([TSl) . There- 
fore, all atoms in the symmetric mode are only involved 
in the dynamics of the atom-photon interactions if the 
system starts with the state |^'i(0)) in Eq. (fT9|) . In fact, 
there are only a few excitations in the asymmetric mode 
due to the atomic interactions. The effect of the exci- 
tations from the asymmetric mode to the dynamics of 
atom-photon interactions is very small. It is because 
the Rabi coupling strength cannot be greatly enhanced 
with a small number of atoms in the asymmetric mode. 
We briefly discuss the validity of this assumption in Ap- 
pendix A. 

It is instructive to express the Hamiltonian in terms of 
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angular momentum operators: 



(+) 



S 



.9+e+, 

e+gl 



The Hamiltonian can be rewritten as 



(20) 
(21) 

(22) 
(23) 



For simplicity, we have assumed that the tunneling 
strengths Jg and Jg are equal. We also have omitted 
the constant term hNA/2. 

By a pply ing the Holstein-Primakoff transformation 
(HPT) [30|, the angular momentum operators can be 
mapped onto the harmonic oscillators which are given 

by 



where j = L, R. Now the Hamiltonian is rewritten as 
H2=hY^ {A + S)S,, + hg{aSj+ + ll.c.) + hxS%, (32) 

j=LM 

where S = {Uee ~ Ugg)N/2 and x = C/ee + U, 



We have omitted the constant term h{Uee t ^gg 



gg 2,Ueg- 



2f/eg)7VVl6 + hNA/2 in Eq 



We consider all atoms at the state \g) are initially pre- 
pared in the ground state of the Hamiltonian in Eq. p2p , 
which can be described by a product of two number states 
as 

\^'2m^\N/2)g,\N/2)g,. (33) 

Without loss of generality, we assume that TV is an even 
number. 

Wc apply the HPT such that the angular momentum 
operators can be mapped onto the harmonic oscillators 



S 
Si 



+ 

(+) 



N 



= b'fb 



(24) 
(25) 
(26) 



In the low degree of excitation, the mean excitation num- 
ber (b^b) are much smaller than the total number of 
atoms N. The angular momentum operators can be ap- 
proximated by the bosonic operators [3,[3l|. The effective 
Hamiltonian can be obtained as 



H^^^ = hAb^b - 



%\/7V(afo^ + H.C.). 



(27) 



Note that the effective Rabi frequency is enhanced by 
a factor of ^/N. This effective Hamiltonian H^^ in 
Eq. (|27p describes the interactions between the collective- 
excitation mode and the single mode of the photon field. 



B. Weak-tunneling regime 

Now we investigate the system in the weak tunnel- 
ing regime, where the atom-atom interaction strengths 
are much larger than the tunneling strengths, i.e. 



Ue,Ug,Ueg^J^ 



Jg. In this limit, we assume that the tun- 



neling between the two condensates is effectively turned 
off. The total Hamiltonian can be approximated as 



i?2 = hA{e\eL + e\eii) + hg[a{e\gL + e\gn) + H.c] 

+hU,,[{eleLf + [e'^ReRf] + Wgg[{g{gL? 

+{9r9r?] + 2ft[/eff(4eLgl5L + e^jeflgjj(7fl,).(28) 

Here we have ignored the terms of the tunnel couplings. 

This Hamiltonian can be expressed in terms of the an- 
gular momentum operators: 



(29) 
(30) 

(31) 



St.+ — 



N/2- 



c^c, 



Slz 



c'c - 



S 



R+ 



d\ N/2 - dU, 



S 



Rz 



Sl- = 
N 

~ T' 

Sr- 
N 



:0V/2~cic, (34) 
(35) 



d\ N/2 



dU, (36) 
(37) 



If the mean numbers of the atomic excitations, (c' c) and 
(d'^d), are much smaller than the number of atoms N/2 
in each well, then the Hamiltonian can be approximated 
dllll as 



(2) _ 
cff ~ 



hAyj{c'c + 

where A^ = 2A + d- 



dU) + hg^j^[a{c) + d^) -f H.c] 

+ {d^dn (38) 
XN /2. The effective Rabi frequency 
is enhanced by a factor of \J N /2. The parameter x is 
much smaller than the effective Rabi frequency because 
the scattering lengths of the inter- and intra-component 
condensates of ^''Rb are very similar [6] . We will ignore 
the terms with the parameter x Eq. (j38|) in our later 
discussion. 



The effective Hamiltonian Hl^J in Eq, 



describes 



the interactions between the single mode of the photon 
field and the two modes of the collective excitations of the 
atoms in the left and right potential wells, respectively. 
This system can be described by a system of three cou- 
pled harmonic oscillators. The effective Rabi frequency 
for each atomic mode is proportional to the factor \J N/2. 
This is different to the effective Rabi frequency, in the 
strong-tunneling regime, which is proportional to the fac- 
tor Jn. 



DARK STATES AND QUANTUM 
DYNAMICS OF THE SYSTEM 



IV. 



We now study dark states of the system which has 
different dark-state subspaces in the strong- and weak- 
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FIG. 2. (Color online) Time evolution of the mean photon 
number (a) and the mean atomic excitations (b) with the 
damping rate k = 100(? and the detuning A = 0. The different 
number of atoms A'^ are 5 x 10^ (black-solid line), 1 x 10"* (blue- 
dashed line) and 2 x 10'* (red-dotted line), respectively. 



tunneling regimes. Let us first introduce the definition 
of dark states. Dark states [2l| are the eigenstates of the 
atom-photon interaction operator V, with zero eigenval- 
ues, i.e., 

V|Dark) = 0|Dark), (39) 

= 0. (40) 

Dark states, in the strong- and weak-tunneling regimes, 
in this system can be obtained as 

h!,^^\D),=0, (41) 



where H^g are 



the two effective Hamiltonians in 
Eqs. dSll) and dSH]) with zero detunings (A = A„ = 0) 
and j = 1, 2. 

In the limit of strong tunnel coupling, the dark state 
\D)i is the product state of the vacuum state of the pho- 
ton field and the ground state of the atomic mode 6, 
which is given by 



\D)i - |0),|0),. 



(42) 



This state is the ground state of the coupled system of 
the atoms and the photon field. 

In the weak-tunneling regime, the system has a family 
of dark states. The family of dark states are 



\Dn)2 - \0)am 



where 



\D?,) 



l)\q^\n-j),\j)a, 



(43) 



(44) 
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FIG. 3. (Color online) Dynamics of the mean photon number 
and mean atomic excitation numbers with the damping rate 
K = lOOg and the detuning A™ = 0. (a) the mean photon 
number (a^a) as a function of the time gt. Time evolution of 
the mean atomic excitation numbers of the atomic mode c, in 
(b), and the atomic mode d, in (c) as shown. The different 
number of atoms A'^ are 5 x 10'^ (black-solid line), 1 x 10* 
(blue-dashed line) and 2 x 10* (red-dotted line), respectively. 



and C" is the binomial coefficient. The dark states |£'„)2 
are the product state of the vacuum state \0)a of the 
photon field and the states ID") are the eigenstates of 
the operator c + d with zero eigenvalues. Note that the 
states in Eq. is a superposition of the states |n— 
j)c\j)d which have the same degree of atomic excitations. 

To gain more insight into dark states, let us first inves- 
tigate the atom-photon dynamics subject to the dissipa- 
tion of the photon field. For a superconducting resonator 
with the frequency ^ 40 GHz can be cooled down to low 
temperatures 25 mK) [32]- This allows us to consider 
the cavity field being weakly coupled to the reservoir at 
the zero temperature [1^ . Note that the relaxation time 
(several ^s) of the single photon inside the superconduct- 
ing resonator is much shorter than the coherence time 
(^ Is) of the cold atoms [l^, [l^l- The effect of the dissi- 
pation of the atoms caused by the noise of the surface of 
the superconductor is negligible [s^] • The main source of 
the dissipation is the damping of the photon field. The 
dynamics of the system can be described by the master 
equation, for the zero temperature, as [s^ [3^ 

Pj = ~j^[H^es\Pj] + -^(SaPj-a^ - a^aPj - p^a^a), (45) 
where pj is the density matrix of the total system, and 
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j = 1,2. Obviously, the dark states \D)i and |D„)2 
are the steady-state solutions of the master equation in 
Eq. (|45| . Thus, the dark states are robust against the 
dissipation of the photon field. In the strong tunnel- 
ing regime, the steady state is the dark state \D)i. In 
the weak tunneling regime, the state of the condensates 
evolves as a mixture of dark states |I?„)2 through the 
dissipation of the photon field. 

Now we study the dynamics of the system in the 
strong-tunneling regime, where the state is prepared as 
|0)a|l)b and |l)f, is a number state. We plot the time of 
evolution of the mean photon number and mean atomic- 
excitation number in Fig. [2] The mean photon number 
and mean atomic excitations undergo a few oscillations 
and then both of them decay to zero. We also see that 
the faster rate of oscillations can be obtained if a larger 
number of atoms N are used. 

We proceed to investigate the atom-photon dynamics 
in the weak-tunneling regime. The system is initially 
prepared as the state |0)a|l)c|0)d, where |l)c is a number 
state. In Fig. [3l we plot the mean photon number, and 
the mean excitation numbers of the two atomic modes 
versus the time. When the atom-photon interactions are 
turned on, the excitation number of the atomic mode 
c decreases while the mean photon number increases as 
shown in Fig. [31 Afterwards, the mean excitation number 
of the atomic mode d starts to increase. This means that 
the energy of the atomic mode c transfers to the photon 
field and the atomic mode d absorbs the energy from 
the photon field. In this way, the two atomic modes 
exchange the energy via the photon field. The faster rate 
of exchanging energy between the atoms and the photon 
field can be attained if a larger number of atoms N are 
used. We also note that the mean photon number in 
Fig. [31[a) is about half of the mean photon number in 
Fig.lHa). It is because the atoms in the atomic mode d, 
in the weak tunneling regime, absorbs the energy from 
the photon field. 

In Fig. [3l^a), the mean photon number decays to zero 
after a period of time. However, the mean excitation 
numbers of modes c and d remain non-zero as shown in 
Figs. 13] (b) and (c). It is because the state of the atoms 
evolves to a mixture of dark states |-Do)2 and |-Di)2, and a 
single excitation is shared by the atoms in the dark state 
|Di)2- This results in the non-zero excitation numbers of 
the two atomic modes. 



V. GENERATION OF ENTANGLEMENT 
BETWEEN TWO SPATIALLY SEPARATED 
CONDENSATES 

We have shown that the system has the different dark- 
state subspaces in the two tunneling limits. Now we 
study the entanglement between the condensates in the 
two different potential wells in the weak tunneling regime. 
In this regime, the system has a family of dark states 
which can be used for generating entanglement. Here 




0.15 



0.05 0.1 0.15 

FIG. 4. (Color online) Time evolution of the entanglement 
witness in (a) and logarithmic negativity in (b), for the damp- 
ing rate k = lOOg and the detuning A^, = 0. The different 
number of atoms A'^ are 5 x 10'^ (black-solid line), 1 x 10* 
(blue-dashed line) and 2 x 10* (red-dotted line), respectively. 

we consider the tunneling between the wells to be effec- 
tively turned off. Therefore, the two independent conden- 
sates in the two potential wells are initially unentangled. 
We will show that steady-state entanglement between the 
two condensates can be produced by evolving to a mix- 
ture of dark states! |Z)„) 2} through the dissipation of the 
photon field |22424l |. 

To study the quantum entanglement between the two 
atomic modes c and d, it is necessary to obtain the den- 
sity matrix of the atomic condensate. By tracing out 
the system of the photon field, we can obtain the density 
matrix pcd, 

Pcd = Tva{p), (46) 

where p is the density matrix of the total system. Let 
us first examine the entanglement of a single dark state 
\Dn)2- For a dark state |D„)2 in Eq. the density 

matrix pcd is given by 

p,rf = |Z?^)(Z?,"J, (47) 

where jZJJJ) is the state in Eq. The degree of entan- 

glement between the two atomic modes can be quantified 
by the von Neumann entropy. It is defined as 

EF^~Tr{pc\npc), (48) 

where pc = Trd,{pcd) is the reduced density matrix. The 
von Neumann entropy is 

n 

Ef = -2-" ^ In |2-"C77|. (49) 
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FIG. 5. (Color online) Plot of the dynamics of entangle- 
ment, (a) entanglement witness W and (b) logarithmic neg- 
ativity Ej^{pcd) as a function of the time gt. The initial 
state |0)a|n)c|0)d with the different excitation numbers n are 
shown, for n = 1 (black-solid line), n = 2 (blue-dashed line) 
and n = 3 (red-dotted line) , respectively. The parameters are 
K = lOOg, A„ = and Af = 5xlO^ 

Thus, the state |Z?JJ) is an entangled state. The degree 
of two- mode entanglement becomes higher for larger n. 

In general, this density matrix pcd is a mixed state. To 
quantify the entanglement of a mixed state, the logarith- 
mic negativity can be used [s^. The definition of the 
logarithmic negativity is [s^ 

i?^(Pcd) -log2 II II, (50) 

where is the partial transpose of the density matrix 
Pcd and II • II is the trace norm. 

However, the logarithmic negativity is difficult to be 
experimentally determined. It is very useful to study an 
experimentally accessible quantity to detect the quantum 
entanglement between the two bosonic modes |37| . If an 
inequality 

|(cdt)|2 > (51) 

is satisfied [37|, then the state is an entangled state. Here 
TT-c = c^c and Ud — d^d are the number operators of the 
atomic modes c and d, respectively. For convenience, this 
quantity W is defined as 

W-Kn<i)-|(cdt)|2. (52) 

If W is negative, then the state is non-separable. This 
quantity W is called as an entanglement witness [38|. 

We investigate the dynamics of entanglement between 
the two atomic modes. We consider an initial state as a 



product state of the three modes, i.e., |0)a|l)c|0)d, where 
|l)c is a number state. We plot the entanglement wit- 
ness and logarithmic negativity versus time as shown in 
Fig. 21 This figure shows that the entanglement wit- 
ness decreases and logarithmic negativity increases with 
a similar rate, and then they saturate after a short time. 
This shows that the steady-state entanglement can be 
produced in a short time via the dissipative photon field. 
The entanglement can also be produced faster if a larger 
number of atoms are used. Besides, we can see that the 
entanglement witness is consistent with the logarithmic 
negativity to indicate the degree of entanglement. The 
entanglement witness is a faithful indicator for detecting 
the entanglement between the two bosonic modes. 

Next, we study the generation of entanglement by us- 
ing an initial state |0)a|?^)c|0)(J with a higher degree of 
excitation, where |n)c is a number state and n is larger 
than one. In Fig. [SJ the entanglement witness and loga- 
rithmic negativity are plotted versus the time. It shows 
that a higher degree of the entanglement can be obtained 
if higher excitation numbers n = 2, 3 are used. 

VI. SUMMARY 

We have studied a two-component condensate in a 
double-well potential, where the atoms are magnetically 
coupled to a single-mode of the photon field inside a su- 
perconducting resonator. The system has the different 
dark-state subspaces in the strong- and weak-tunneling 
regimes, respectively, and it gives rises to the different dy- 
namics of atomic excitations in the two regimes. Steady- 
state entanglement between the two spatially separated 
condensates can be produced by evolving to a mixture of 
dark states through the dissipative photon field. We have 
shown that the entanglement can be faithfully indicated 
by an entanglement witness. 
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Appendix A: Validity of the effective Hamiltonian in 
the strong-tunneling regime 

In this appendix, we examine the validity of the effec- 
tive Hamiltonian in Eq. (|27p in the limit of strong 
tunnel coupling. We express the Hamiltonian in term of 



8 




FIG. 6. (Color online) Level scheme of the atoms in the 
double-well potential. In the strong-tunneling regime, the 
atoms, with the two states |e+) and \g+), in the symmetric 
mode are coupled to the cavity field. The two ground states 
of the symmetric and asymmetric modes {\g+) and \g-)) are 
coupled to each other via the atom-atom interactions. 



the symmetric- mode and asymmetric-mode operators as: 

-9-9-) + f^9{ae^+9+ + H.c.) + hg{aeLg- + H.c.) 

^2l 



+ ^ [{gW + 9l9-f + + 9l9+n 

+hUeg [(e^e+ + ele-)iglg+ + gig-) + (e^ 
+ele+)iglg^+gLg+)]. 
Let us define 

F+=glg+, F-=glg_, 



-,{9-9- -9+9+)- 



(Al) 

(A2) 
(A3) 



The commutation relations [F-i,F±] — ±F± and 
— 2Fz are satisfied. Tfie operators F± and F^, 
and S*!"^^ and in Eqs. ^ to 1^ generate a SU(3) 
algebra [s^. In the limit of large N , we can apply the 
HPT to the operators: 



F+ ^f^^N-pf, F- = f^N-pf, (A4) 
Fs = Pf - N/2. (A5) 

Assume that the mean excitation number (/^/) is much 
smaller than N, we can approximate the operators as 
[39f: 



1 



N 



f 



1 



(A6) 



In the low-degree-of-excitation regime, the approximated 
Hamiltonian can be written as 

H « huaa^'a + hw'b% + hgVN{ab'' + H.c.) + H' . 

(A7) 

The Hamiltonian H', contains the terms of the operators 
in the asymmetric mode and the terms from the nonlinear 




FIG. 7. (Color online) The expectation value {f\f) as a 
function of the time Jgt. Different strengths of UggN are 
shown: UggN = Jg (red-solid line), UggN = 5Jg (blue-dashed 
line), and UggN — 10 Jg (black-dotted line). 



interactions, and the constant terms are omitted, which 
can be written as 



H' 



-eLe_) -|- hg{ae^_g- -f H.c.) 



Ji(A4 



./e)(eLe_ 
■[{eW 



+e^^e-f + {e\e- + eLe+)2] + 
+9-9-f + f^Ueg [{ele+ + ele^){glg+ + glg^) 



N{e\t 



eLe+)(/t + /)]. 



(A8) 



Here we consider the number of atoms in the excited 
states to be very small. We also assume that the strength 
of the Rabi coupling g is weak compared to the tunneling 
strength Jg and nonlinear strength UggN but g is much 
stronger than Uee,Ugg and Ueg. Therefore, the Hamil- 
tonian H' can be approximated by the Hamiltonian H" 
as 



where 



Xl=Jg + UggN, 

\2 



UggN 



(A9) 

(AlO) 
(All) 



From Eq. (|A9p . nonlinear interactions can give rise to 
the transitions of the atoms in the symmetric mode to the 
atoms in the asymmetric mode, and vice versa. The level 
scheme is shown in Fig. [51 Note that this Hamiltonian 
H" is exactly solvable. The time-evolution operator can 
be factorized as 1351 



Pt2 



(A12) 



= exp(-f 



exp 



In(Ai) 



exp 



A 



(A13) 
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where 



vacuum state |0)/ of the mode /. The state becomes 

(A18) 



i,,(i))=Arf;4fn(i^)»i2„),. 



Ai = ( cosh/3 
A2 



n=0 



2/3 
2A2 sinh /3 



sinh j3 



2/3 cosh /3- A'l smh/3' 
-2iXit, A2 = -2iX2t. 



(A14) 
(A15) 

(A16) 
(A17) 



The mean excitation number (f^f) is 
(/V) = lA 



i/4|2^ n(2n) !A^" 

2^ 22ji-l 



2 

(n!)2- 



(A19) 



In Fig. [71 we plot the expectation value (/V) versus 
the time, for the different strengths of atomic interaction 
UggN . It is shown that there are only a few excitations 
in the asymmetric mode even if the atomic-interaction 
strength UggN is much larger than the tunneling strength 
Jg. Therefore, the Rabi coupling strength cannot be 
We then apply the time-evolution operator S{t) to the greatly increased due to the collective enhancement. 
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